Multiscale Mixing Efficiencies for Steady Sources 



Charles R. Doering* 

Department of Mathematics and Michigan Center for Theoretical Physics, 
University of Michigan, Ann Arbor, MI 48109-1043, USA 



o 
o 

(N 

< 



i 

O 

Oh: 



(N : 

> ■ 

^ ■ 

in 
oo 

o 
o 

^ : 

Oh. 



X 



Jean-Luc ThifFeault^ 
Department of Mathematics, Imperial College London, SW7 2AZ, United Kingdom 

Multiscale mixing efRciencies for passive scalar advection are defined in terms of the suppression of 
variance weighted at various length scales. We consider scalars maintained by temporally steady but 
spatially inhomogeneous sources, stirred by statistically homogeneous and isotropic incompressible 
flows including fully developed turbulence. The mixing efficiencies are rigorously bounded in terms 
of the Peclet number and specific quantitative features of the source. Scaling exponents for the 
bounds at high Peclet number depend on the spectrum of length scales in the source, indicating 
that molecular diffusion plays a more important quantitative role than that implied by classical 
eddy diffusion theories. 
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Introduction. Mixing in fluid flows plays a central role 
in many scientific and engineering applications and is the 
subject of a large body of theoretical research [1]. In this 
paper we consider the mixing of a passive scalar main- 
tained by a temporally steady but spatially inhomoge- 
neous source without a clear separation of length scales 
between the source and the stirrer, as is necessary for 
multiscale analyses, or when length scales in the source 
are the smallest scales in the problem. 

The effectiveness of a given stirring velocity vector field 
in a statistically stationary state is naturally gauged by 
its ability to suppress the space-time averaged variance 
of the scalar: a decrease in the variance indicates a more 
uniformly mixed scalar. The ratio of the scalar variance 
in the presence of stirring to that resulting from molec- 
ular diffusion alone defines a dimensionless measure of 
the mixing efficiency of the flow. Scalar variances may 
be weighted at diverse spatial length scales, so we in- 
troduce a family of mixing efficiencies that measure the 
effectiveness of the stirring on different scales. Not un- 
expectedly, these mixing efRciencies generally depend on 
details of both the source and the stirring. Sawford & 
Hunt [2] highlighted the role of source size in particle 
dispersion models for the case of small sources, showing 
that the concentration variance depends explictly on the 
molecular diffusivity. Here we make direct use of the 
advection-diffusion equation to show the nature of this 
dependence for a wide range of source types. 

For statistically stationary, homogeneous and isotropic 
flows — including turbulence — we derive new, mathemat- 
ically rigorous and physically relevant limits on these ef- 



ficiencies in terms of the Peclet number and quantitative 
features of the source alone. We discover that the high- 
Peclet number scaling exponents of the efficiency bounds 
depend on the spatial dimension, the length scales upon 
which the variance is observed, and spatial scales in the 
source. These results show that the spectrum of length 
scales in the source may be more important than the 
distribution of length scales in the flow in determining 
the mixing efficiencies. Moreover, we derive novel quan- 
titative estimates of nontrivial scaling exponents for the 
mixing efRciencies in cases where the source is fractal or 
measure-valued. These are anomalous scalings that can- 
not be deduced from dimensional arguments alone. 

Problem statement and definitions. The advection- 
diffusion equation for the concentration 6(x,t) of a pas- 
sive scalar is 



(1) 



where k is the molecular diffusivity. We restrict attention 
to spatially periodic boundary conditions for x e [0, L]''. 
The time-independent source s{x) in (1) is taken (with- 
out loss of generality) to have spatial mean zero so 
that eventually the concentration 9{x, t) will have spa- 
tial mean zero as well. 

The velocity field u{x, t) is given. It could be the 
solution of some dynamical equations or a particular 
stochastic process, but in any case we consider it to be 
a prescribed divergence-free vector field with the follow- 
ing equal-time single-point statistical properties shared 
by homogeneous isotropic turbulence: 
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where overbar represents the long-time average (assumed 
to exist) at each point in space. The r.m.s. velocity U 
measures the strength of the stirring and ^l indicates the 
flow field's strain or shear content. The ratio X = U/Q 
corresponds to the Taylor microscale for homogeneous 
isotropic turbulence. The Peclet number for the flow is 
Pe = UL/k. 

We quantify the mixing of the scalar by the magni- 
tude of the variances (|V6'p), (6*2), and (|V~^0p) where 



a bound that depends on the "shape" of the source func- 
tion but not its amplitude, and on the stirring veloc- 
ity field only through its influence on the length scale 
k/U = Pe-^L. 

Limits for the small scale and large scale efficiencies 
£±i are obtained from (4) in the same manner after inte- 
grations by parts and application of the Cauchy-Schwarz 
inequality. For the gradient variance 



(•) denotes space-time averaging. (The operator V'^ is i^Ps) = {{uy^ + «V^) • V9) < {\uip + KVip\^) {\V9\^) 



defined by its Fourier representation in the periodic do- 
main, — ife/fc^.) These variances measure the 
fluctuations of 9 at relatively small, intermediate and 
large length scales. Collectively we write the norms 
(|V2'6'|2) for p = 1, 0, and -1. Note that {\W-^e\^) is 
reminiscent of the recently-introduced "mix- norm" [3], 
as both downplay the importance of the small scales. 

In order to define dimensionless mixing efficiencies 
we use the baseline variances defined by the solution 
9q of equation (1) with the same source but u = 0, 
i.e., ^0 = — K~^A~^s. Comparing fluctuations in the 
presence of stirring to the moments (jV^^oP) allows us to 
gauge the effect of stirring. We deflne the dimensionless 
mixing efficiencies £p by 



(|VP^?oP)/(|vf^^|2) 



(3) 



These efficiencies increase when the stirring decreases the 
scalar variances relative to molecular diffusion alone. 

Mixing efficiency bounds. Upper bounds on £p result 
from lower limits on the variances. We begin with esti- 
mates on (^^^) using the method developed in [4]: multi- 
ply (1) by an arbitrary (but smooth, spatially periodic) 
function (f{x), average, and integrate by parts to find 

{9{u-W + kA)<p) = - (ifis) . (4) 

Fluctuations are bounded from below via the min-max 
variational expression 

((9^) > maxmin{(^?2) | {u ■ V <fi + k Aip)) = - {<ps)} . 

The minimization over "& is easily achieved, equivalent to 
application of the Cauchy-Schwarz inequality, yielding 

(61^) > max (ipsf /{{u-Wip + KAipf) . (5) 

This is the sort of variational estimate derived in [4]. 
Plasting and Young recently enhanced that analysis by 
including the scalar dissipation rate as a constraint [5]. 

The subsequent maximization over ip is particularly 
simple for statistically homogeneous and isotropic flows 
satisfying (2), for then the denominator in (5) is 

{{u ■ + «A^)2) = (^|V<^P + «2 (A<^)^^ , (6) 

i.e., the quadratic form ((^[k^A^ — {U^ / d) A](fi) . Hence 
the variational problem (5) yields 



(sA-^s) 



{s{A^-{UyK^d) A}-^s) ' 



(7) 



so 



{\V9\'^) > max {ifsf / {{uip + nV^f) . (8) 

A potentially sharper bound involving the full two-point 
correlation function for the velocity field can be ob- 
tained by formally minimizing over 9 [6], but for our 
purposes the estimate (8) suffices. For statistically 

homogeneous isotropic flows the denominator above is 
(<^[— K^A -|- U'^](f) and optimization over (p leads to 



(,s(-A-» 



s{-A + U^/k^}-^ s 



(9) 



The bound on large scale fluctuations follows from (4) 
using 9 = V- V~^^, an integration by parts, and Cauchy- 
Schwarz: 

{ipsf = (V(m -Vip + KAip)- {V-^9) f 

< {\Vu -Wip + u- VV<p + kAV<p\^) <|V-i6'p) 



so that 



max 



{psY 



f (|Vu • V(p -(- u • VV^J -h kAV(p|2) 



For statistically homogeneous isotropic flows the denom- 
inator is ((^[-K^A^ (fJ'^ld) A^ - {n'^/d) A\>p>) so that 



{s{-A-^)s) 



{-A3 + {V^/K^d) A2 - {n-^/K^d) A}-^ s 
It is helpful to rewrite the bounds in Fourier space: 



£\< 



Efel«feW + Pe')' 

Efc|sA=IV(fc^ + Pe'A:Vd)' 

EfckfelVfc'^ 

EfekfelV(fc^ + Pe^A;V<^ + Pe^fcVA2d) 



(10a) 
(10b) 
, (10c) 



where we have rescaled [0, L]'' to [0, 1]'' so that wavevec- 
tor components are integer multiples of 27r. Now we 
investigate the large Peclet number behavior of these 
bounds for a variety of classes of sources. 
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Monochromatic Sources. For sources that depend only 
on a single wavenumber k^, the bounds simplify to 



fi < ^l + PeVfcf, 
So < ^l + PeVfc2d, 

£_i < ^J 1 + PeVfc^d + Pe'^/X^k^d . 

Each efficiency is asymptotically proportional to Pe, cor- 
responding to the expected suppression of variance if the 
molecular diffusivity k is rc^placcd by an eddy diffusivity 
proportional to UL. Moreover these upper bounds are 
sharp: they may be realized by uniform flow fields whose 
direction varies appropriately in time to satisfy the weak 
statistical homogeneity and isotropy conditions used in 
the analysis [5, 6] . Each estimate also exhibits a decreas- 
ing dependence on the length scale of the source: at high 
Pe the estimates for the small and intermediate scale ef- 
ficiencies £i and So are ~ Pe/fcg. This suggests that an 
eddy diffusivity might better be defined as a product of 
U and a length scale ^ k~^ characterizing the source, 
rather than some scale characterizing the stirring. 

Square Integrable Sources and Sinks. In situations 
where the Fourier coefficients are square-summable, the 
asymptotic Pc ^ oo behavior of the mixing efficiency 
bounds for smooth sources is straightforward: 



Si<Pe. 



SQ<Pe, 



I EM' 



S-i < Pe . 



(lla) 
(lib) 
(11c) 



These are the same Pe scalings as observed for 
monochromatic sources but with prefactors involving dis- 
tinct characteristic length scales of the source tailorcid to 
the scales where the different efficiencies are tuned. Wc 
reiterate that the Pe^ scaling of the efficiencies is precisely 
that which is expected from the conventional notion of 
eddy diffusion, at least with regard to the U and k depen- 
dence. A novel feature of these rigorous estimates is that 
the "mixing lengths" to be employed in constructing the 
effective diffusion from U depend on (i) the source struc- 
ture rather than some characteristic persistence length in 
the flow, and (ii) the length scales in the concentration 
fluctuations stressed by the different multiscale efficien- 
cies. 

Rough Sources. The Peclet number scaling may actu- 
ally change for "rough" sources, i.e., when s{x) is not 
square integrable, resulting in anomalous behavior for 
some of the efficiencies. The roughest physically mean- 
ingful sources are measure- valued sources like (5-functions 
with nondecaying Fourier coefficients \sk\ = 0{1) as k 
oo. Then the sums in (lla) and the denominator of (lib) 



diverge in = 2 or 3 so those scalings are invalid. In this 
case the Pe dependence of Si disappears altogether and 
all finite kinetic energy stirring fields are completely in- 
effective at suppressing small scale fluctuations. 

To determine the high-Pe behavior of So we approxi- 
mate sums by integrals. The denominator of (10b) is 



k'^-'^ dk 



k^ + {Peyd)k'^ J2^ k^ + Pe^k'^/i^n'^d) 



For d = 2 the integral in (12) is 
fcdfc 



/ 

J 271 



log Pe 



fc4 + Pe2fc2/87r2 Pe^ 



resulting in the asymptotic bound 



^o<Pe/VlogPe, d = 2. 



(12) 



(13) 



(14) 



Hence there is at the very least a logarithmic correction 
to So as compared to the square integrable source case. 
For d = 3 the integral in (12) becomes 



fc^ dfc J_ 

V+p^Fkyu^ ^ p^ 

resulting in an anomalous scaling bound 



So < \^ , d = 3. 



(15) 



(16) 



This is a dramatic modification of the classical scaling. 

Similar analysis shows that the upper bound on the 
large scale mixing efficiency S-i ~ Pe in (11c) persists 
even for these roughest sources. 

Rough sources with a cutoff. Approximate delta-like 
sources of small but finite size i have Fourier coeffi- 
cients Sk that arc approximately constant in magnitude 
up to a cutoff wavenumber of order 27r/£, beyond which 
the spectrum decays as for a smooth source. We may 
deduce the behavior of the bound on So for such sources 
by inserting an upper limit at L/£ ^ I into the integral 
in (12). For large but intermediate Peclet numbers sat- 
isfying 1 <C Pe <C L/£, the cutoff is ineffective so the 
logarithmic correction (14) in d = 2 and the anomalous 
scaling (16) in ci = 3 appear. However for Pe ^ L/£, i.e., 
when the modified Peclet number based on the smallest 
scale in the source UI/k^ 1, the smooth source results 
apply and we recover the mixing efficiency bounds linear 
in Pe, as in (11). Figure (1) shows this scaling transition 
for the d = 3 case. Even in the ultimate regime where 
the source appears smooth, the prefactor in front of the 
high-Pe scaling bounds are significantly diminished by 
the small scales in the source: So ^ [log{L/ £)]~^^' Pe in 
d 2, and So < [e/L]^/^Pc in = 3. 

Fractal Sources. We may also analyze anomalous scal- 
ings for more general "fractal" rough sources where the 
Fourier spectrum |sfc| decays as k~'^ with < 7 < d/2. 
The roughest measure-valued sources have 7 = while 
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FIG. 1: Upper bound for the mixing efficiency £o as a func- 
tion of Peclet number for a small source with £ = 10~^L 
stirred by a three-dimensional statistically homogeneous and 

isotropic flow [computed from Eq. (10b)]. The intermediate 
Pe^^^ scaling for 1 <C Pc <C (L/£) is evident. 

for 7 > d/2 the source is square integrable and thus effec- 
tively smooth as far as these niultiscale mixing efficien- 
cies are concerned. In order to examine the high Peclet 
number asymptotics of the bounds on the various £p we 
estimate integrals similar to (12) but with an extra fac- 
tor of A:"^"*" in the numerator arising from |sfep. The 
results are summarized in Table I. In d = 2 the scal- 
ing for £i is anomalous for any degree of roughness while 
£o is anomalous only for 7 = 0. In d = 3, f 1 is again 
anomalous for any degree of roughness while £0 scales 
anomalously for < 7 < 1/2. For both d = 2 and 3 
the bound on the large scale mixing efficiency £-1 is al- 
ways conventional. Of course these scalings neglect any 
large-A; cutoff for the rough sources. If there is a cutoff 
at wavenumber 27r/^ then the same arguments apply to 
recover the normal scaling ^ Pc for Pc^ L/£. 

Discussion. The multiscale mixing efficiency bounds 
derived here reveal new aspects of the ability — and in 
some cases the inability — of statistically homogeneous 
and isotropic incompressible flows to effectively suppress 
fluctuations in passive scalars. Several notable features 
of the phenomena have emerged from this analysis. One 
is that the structure of the source maintaining the scalar 
concentration plays a central role in efficiency of mixing 
at different length scales while the detailed structure of 



the flow, i.e., the spectnun of length (and time) scales 
in the velocity field, plays only a secondary role as far as 
the small and intermediate scale efficiencies £1 and £0 are 
concerned. The bound on f_i in (11c), however, means 
that small scale structure and strain in the advecting flow 
could possibly enhance the large scale mixing efficiency. 

For rough sources the efficiencies scale at least as 
anomalously as indicated in Table 1, and no clever small 
scale stirring can alleviate this effect. When efficiency 
bounds scale anomalously, molecular diffusion plays a 
much more important role than implied by conventional 
eddy diffusion theories. For then — as far as variance sup- 



TABLE I: Scalings of the bound on the mixing efficiency £p 
as functions of the source roughness exponent 7 of the source 
in two and three dimensions. 
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Pe 
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Pe 
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Pe 
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pression is concerned there can be then no "residual" 
effective diffusivity due to stirring in the limit of negligi- 
ble molecular diffusion. 
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